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HOMOGENIZATION OF BRINKMAN FLOWS IN HETEROGENEOUS 

DYNAMIC MEDIA 
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Abstract. In this paper, we study Brinkman’s equations with microscale properties that are 
highly heterogeneous in space and time. The time variations are controlled by a stochastic 
particle dynamics described by an SDE. The particle dynamics can be thought as particle 
deposition that often occurs in filter problems. Our main results include the derivation of 
macroscale equations and showing that the macroscale equations are deterministic. The 
latter is important for our (also many other) applications as it greatly simplifies the macroscale 
equations. We use the asymptotic properties of the SDE and the periodicity of the Brinkman’s 
coefficient in the space variable to prove the convergence result. The SDE has a unique 
invariant measure that is ergodic and strongly mixing. The macro scale equations are derived 
through an averaging principle of the slow motion (fluid velocity) with respect to the fast 
motion (particle dynamics) and also by averaging the Brinkman’s coefficient with respect to 
the space variable. Our results can be extended to more general nonlinear diffusion equations 
with heterogeneous coefficients. 
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1. Introduction and formulation of the problem 

1.1. Motivation. In many porous media application problems, the media is subject to change 
due to pore-scale processes. For example, in filter applications da m mi [n], the media 
properties and microscale geometry change due to particles are captured by the filter (see 
Figure for illustration). In this figure, we depict a filter element and particle deposition 
process (following [IS])- The particle deposition changes the microscale geometry of the filter 
and thus can greatly affect its macroscopic properties that are used in simulations dais]- 
The change due to particle deposition is described by Stochastic Differential Equations (SDEs) 
m where the particles’ mean velocities are affected by the fluid velocity. Thus, the modified 
effective properties strongly depend on particle dynamics and deriving and understanding these 
effective properties are essential for many of these applications. Motivated by this application, 
we consider a Brinkman model (cf. |15L IT6] ) where the permeability changes due to particle 
dynamics that are driven by the fluid flow. 

In the paper, we derive macroscopic model assuming that the particle dynamics at microscale 
occurs in a much faster time scale compared to the flow. This is typical in these applications 
due to the particle dynamics and their interaction. We derive a macroscopic model where the 
new upscaled permeability is computed using spatial microscale variations of the permeability 
and fast dynamics. Besides computing the permeability value, we show that the permeability 
is deterministic. This is a useful findings as it allows to compute the upscaled permeability in 
a deterministic manner and avoid stochastic macroscale PDEs. 
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Filter element Filtering medium 



Figure 1. Illustration of a filter element (top) and a cross section of a filter 
with a particle deposition (bottom). Courtesy of G. Printsypar. 


Even though our application is specific for Brinkman’ equations, our mathematical concepts 
can be used for many important applications where the media properties at the microscale are 
affected by SDE. An example can be a diffusion equation with heterogeneous coefficients that 
depend on a field described by SDE. In general, we have 


dv^{t, x) 


1 

e 




(Jl!^ T 

— {t,x) = L{a{-,v%un 
u^{t, x))dt + ^J~^dW (t, x). 


( 1 . 1 ) 


where W{t) is a standard Brownian motion. Eor example, L(a(|, u^), n^) = V • 

The question of interest is to derive macroscale equations which will be investigated in our 
future works. 


1.2. Mathematical model. We consider the following system: 

' / X \ 

— (t, x) = Au^(t, x) - x) - a x)j u^{t, x) + f(t, x) 

div u‘^{t,x) =0 


dv^{t, x) = — x) — u^(t, x))dt + y—dW{t, x) 

u^{t, x) =0 
u^{0,x) = Uq{x) 

u^(0,x) = Vq(x) 


in [0, T] X D, 
in [0, T] X D, 

in [0, T] X D, 

on [0, T] X dD, 
in D, 
in D, 


( 1 . 2 ) 
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where D is a bounded domain of with a smooth boundary dD, and are the velocity 
and the pressure of the fluid and the velocity of a particle, a is the Brinkman coefficient and 
describes the porosity of the medium and is affected by the particle deposition as explained 
in the previous section. W{t) is an L^(Ii))^-valued standard Brownian motion defined on a 
complete probability basis (17, Jr, IP) with expectation E, and Q is a bounded linear operator 
on LP‘{DY of trace class. Uq and Vq are the initial velocities for the fluid and for the particle, 
and / is an external force. 

The first equation in system (1.2) is a Brinkman type equation. For a fixed e and when the 
equation for u‘" is not coupled with the equation for v^, number of papers have been devoted 
to its mathematical analysis see e.g. |18j and the references therein. In [18] the authors 


study the existence, uniqueness and regularity of solutions for the Brinkman-Forchheimer 
equation in dimension 3 and their global attractors. Let us observe that in [T8|, the term 
a{u)u is assumed to be monotone, while it’s not in the current paper. In particular, the 
authors in [T8| studied the polynomial case thoroughly. Their result was extended to the 
convective Brinkman-Forchheimer equation which includes the nonlinear term of the Navier- 
Stokes equations. Similar results are obtained in m where a dissipative term is added to the 
equation. 

Our main goal in this paper is to study the asymptotic behavior of the solutions of system 
Notice that is random through the function a that depends on the 


(1.2) when e 


0 . 


stochastic process solution of a stochastic differential equation. Moreover, the function a is 
a multi scale function. Here, is the slow component and is the fast one. We will prove 
that converges to an averaged velocity u solution of the averaged equation (5.3) where the 
averaged operator a is given by (5.2). Here, the averages are taken with respect to the periodic 
variable y and the invariant measure associated to the process for a frozen u^. 

There is a quite large number of papers dealing with averaging principles for finite dimen¬ 
sional systems in both deterministic and stochastic systems. Less has been done in the infinite 
dimensional setting, we refer to lai and the references therein. There is not much in the 
literature dealing with averaging systems in porous media. In [8|, the authors prove an aver¬ 
aging principle for a very general class of stochastic PDFs. Our system looks similar to theirs 
with a very important difference. Our function a{-,v)u is not Lipchitz and it contains the 
variable x/s that describes the heterogeneities of the medium. Hence, their results although 
very general can’t be used here. 

For e > 0 hxed, the well posedness of system (1.2) does not follow from classical results and 
has to be studied accordingly. In this paper, we assume that a £ C^. Despite this condition, 
the term a{-,v)u in system (1.2) is neither Lipschitz nor autonomous. 

We prove the existence of weak solutions by using a Galerkin approximation (u^, that is 
solution of a well posed system and then pass to the limit on n after performing some uniform 
estimates in n. These estimates are also uniform in e. By using our assumption on a and the 
special form of our system, we are able to prove the uniqueness of the weak solution 
We prove that our weak solution is also strong, and get better uniform estimates in e for the 
solution in the Sobolev space 

Then, we study the asymptotic behavior of the fast motion variable for a frozen slow 
motion variable u^. Indeed, we consider the SDE (4.1) for a given It has a mild solution 
which is also a strong solution. Its transition semigroup is well dehned and has a unique 
invariant measure which is ergodic and strongly mixing. We will consider the auxiliary 
equation 


c(e)T^(r/,^)= / [a^r]) - {^)) . 

Jd 
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for c(e) > 0 to be chosen later. The solution is given by 


roo r r 

Jo Ud 


dt. 


Here is the generator of the semigroup P^, ^,r] ^ L^{D)^ and cj) is an arbitrary function 
in V. The operator is defined in section 3 while the operator is defined in section 5 and 
refers to the averaged of wrt to the invariant measure The main difficulty stands in 
passing to the limit on the term 


L 


{a^{t))u^{t) — a{u{t))u{t)) (pdx 


(1.3) 


for (/> G V, where a is defined in (5.2) at the beginning of section 5. This is done by using 


the Ito formula on and isolating the term (1.3). We mainly follow the idea already 

introduced in |8]. By using the uniform estimates obtained in Section 3, a tightness argument 
and some known results for periodic functions, see [T] (lemma 1.3) the passage to the limit is 
performed in distribution. We obtain a convergence in probability by using the fact that the 
limit u is deterministic. 

The paper is organized as follows, Section 2 is dedicated to the introduction of the functional 
setting and assumptions. In section 3, system (1.2) is analyzed for every e > 0. In particular 


existence of strong solutions are established with their uniqueness and their uniform estimates 
with respect to e. The fast motion is analyzed in section 4, these are known results and 
we give some references. The passage to the limit is performed in section 5. Furthermore, the 
well posedness of the averaged equation is established. 


2. Preliminaries and Assumptions 


We make the following notations for the spaces that will be used throughout the paper. For 
any Hilbert space X, B},{X) denotes the Banach space of bounded Borel functions </> : A —)• M 
endowed with the supremum norm; 

U\\b,{X) = sup |(/.(x)|. 

x&X 

Cf,{X) denotes the subspace of bounded and continuous functions and C^{X) the subspace of 
functions that are k times Frechet differentiable with continuous and bounded derivatives up 
to order k endowed with the norm: 


U\\cpx)= “ax sup 
b'- > 0<l<k^^x 


dx^ 


Ci{X) 


where C^{X) =R and for every i > 0, T®(A) = £(A; £®“^(A)), the Banach space of bounded 
linear operators from X to 0~^{X). For Y = [0,1]^ the space C^{Y) denotes the space of 
continuous functions on Y that are T-periodic and the space L'^iY) denotes the closure of 
C#{Y) in L\Y). 

We denote by H and V respectively the closures of V in L?‘{D)^ and where 

V := |n G C°° I div u = 0, u = 0 on 9Zl| . (2.1) 

H is a Hilbert space with inner product inherited from Lp‘{D)^ denoted by (•,•). Denoting 
by H' and V' the dual spaces, if we identify H with H' then we have the Gelfand triple 
V C H C V' with continuous injections. The dual pairing between V and V' will be also 
denoted by (•,•). The norm in any space X will be denoted by || • ||x 
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The function a is positive and satisfies the following conditions: 

for every a: G 

for every y , with 

\\^{y^ OllcihRS) — 


( 2 . 2 ) 

(2.3) 

(2.4) 

independent of y G T. We remark that condition (2.4) implies uniform boundedness for a on 
Y X as well as Lipschitz condition in the second variable, uniform with respect to the first 
one. 

We also assume that / G L‘^{0,T; L‘^{D)^). 

3. Study of the system 


In this section we prove the existence and uniqueness of the solution for the system (1.2) as 
well as some uniform estimates. 


3.1. Well-posedness of the system (1.2). For any e > 0 we denote by the operator, 

■. {D)^ ^ {D), a^{r]){x) = a . (3.1) 

Let us show that is a well defined operator. Given the condition (2.3), we need only to 
show the measurability in x of a‘^{r]) for any rj G For such a function, we consider 

a sequence rjn G Co{D) convergent to y pointwise in D. The function (y,x) —>■ a{y,r]n{x)) 

is a Caratheodory function, measurable in y and continuous in x, so x —>■ a is 

measurable, and by the Lipschitz condition of a is pointwise convergent to a^{rj), which shows 
that a^(r]) is measurable. Moreover 

Theorem 3.1. Assume that G H for every e > 0, then for each T > 0, there ex¬ 
ists a unique solution of the system (1.2), u'^ G L°°(ll; C([0, T]; H) n L^(0, T; V)) and G 
L^(ll; (^([O, T]; L^(Zl)^) in the following sense: P a.s. 


ID 


'{t)4>dx— f UQfdx + f f Vu^{s)'V4>dxds + ( f {v^)u^fidxds 

JD Jo Jd Jo Jd 


lo Jd 



for every t G [0, T] and every (/> G V, and 

v^{t) = uge-*/" + - f u^{s)e-J-^^/^ds + ^ f e-^^-^'^/^dW{s). 


^ Jo 


Jo 


f{s)(j)dxds, 

(3.2) 

(3.3) 


Moreover, if the initial conditions Uq are uniformly bounded in H, then the solutions 
satisfies the estimates: 


sup ||m^||loo(q.j;^2(o,T;V)) < C't, 
£>0 

sup ||u^||L°om.CQ 0 ^'r].H)) < Ct, 

e>0 


and 


sup 

£>0 




dt 


< Ct- 


(3.4) 

(3.5) 

(3.6) 


L“(O;L2(0,r;V')) 

Also, if the initial conditions Vq are uniformly hounded in Lf [Q.-, Lf we also have the 

estimate for : 

supE sup \\v%t)\\l2(u)3 < Ct- 
£>0 4e[0,T] 


(3.7) 
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Proof. We prove the existence of solutions through a Galerkin approximation procedure. 
We consider {ek)k>o a sequence of linearly independent elements in V n L°°{D)^ such that 
span{ek | /c > 0} is dense in V. We define the n-dimensional space for every n > 0 as 
span{ek \ 1 < k < n} and we denote by n„ the projection operator from V onto V„. 

Let us denote by the following process 

w%t) = ^ r G C7([0, T\,L^{Df). (3.8) 

Jo 

Moreover, there exists a random constant C{ui) almost surely finite such that 

sup ||u;=(t)||^ 2 (m 3 < C'(w). 
te[o,r] 

Now, in order to prove the existence of solutions, we will proceed using a path wise argument; 
we fix cj G n and define the Galerkin approximation 

«(o;),4(o;)) GC([0,r];H„)nL2(0,r;V„)) xC([0,r];n„L2(D)3) 

solution of the following system 


f —^{t)(j)dx+ f V{t)V(pdx + f + w^{t))u^{t)<l)dx = f f{t)<j)dx, (3.9) 

Jd Jd Jd Jd 

for every f G V^, (0) = 

4(0) = 0. (3-10) 

where 

(3.11) 

Then, we pass to the limit on (u®, z^) when n —)• oo. 

We write u%^{uj,t,x) = Yl'l=ial^{uj,t)ek{x) and z^{u},t,x) = Ylk=iK(^^'^)^k{x), and get the 
following system for the coefficients a| and 


r n 


da^ f II. 

/ SkBjdx + y^ CLki^^i) / VckVejdx + 

k=i k=i^ 

y~^ / ak{oj,t)a‘^ [ w{u!,t)+ y^bkiu),t)ek] ekOjdx = / f{t)ejdx, 
k=iJ d V k=i ) Jd 


9bkr 


= -- (^fc - 4) > 'i-<k<n 


= / Uo^kdx, 1 < k < n 

Jd 


= 0, 1 < A: < n 


a|(w,0) 

I 

for each 1 < j < n. We make the following notations: 

ttij = / ei{x)ej{x)dx, hj = / Vei{x)'Vej{x)dx, fj{s)= / f{s,x)ej{x)dx, 

Jd Jd Jd 


(3.12) 


{F^)ij{(jj,t,bi,...bn) = j W ^u)(a;,t) + 


CiCjdx 


and 
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and the system is written with these notations as: 


da 


+ Y ^khj + Y bl, 6^) = 


k=l 

dt 

«|(0) 

m 


k=l 


k=l 


--(K-o-k)^ 1 < /c <n 


= / u^Gkdx, 1 < k < n 

Jd 


= 0, l<k<n 


(3.13) 

for each 1 < j < n. Given the linearly independence of the sequence {ek)k>o, the definition of 
the functions {F^)ij and the Lipschitz condition satisfied by a, the system has unique solution 
(bl) i<k<n in ^^[OjT] for every T > 0. This means that and z^ = v^ — is a 
solution for: 

r f f‘ 

/.Kit I »—r tr / . \ »—r 11 I c / c / . \ c / . \ \ I n / . \ ! 1 


ID 


dt 


dz^ 




<(o) 

= 

. 4(0) 

= 0, 


4 - o > 


(3.14) 


for every (j) £ V„. We take 4> = uf^ in (3.14), use the positivity of a to derive that 
d 


so 


q^\\''Ai\\l'^{D)^ — II/(^)IIl 2(D)3 + ll^nllL2(D)3 

ll^nllL2p)3 < e* (||/||l2(o,T;L2(d)3) + ||Uq||l2(£,)3) , 


snp ||lfnllL“(0,T;L2(£))3) < (7^(1 + HiIqIIh) 
n>0 


(3.15) 


We also obtain that 

pT f 1 

J l|V<||i2(B)nxuds +-|K(r)||22(^^3 < J y^/(t)u^dxdt + -||u^|||2(^)3 ^ 

r e 2 e 1 £ 

Il^^nlli2(£,)nxn(is < 7'||/||L2(0_T;L2(D)3)||u^||i2('£))3 + -||WoIIl2(D)3) 


so 


sup IKI|i,oo( 0 ,r;V) < C'r(l + I^oIIh)- 

n>0 


(3.16) 


The estimates (|3.15|) and (|3.16|) imply using the first equation of the system (3.14) that 

< C't(1 + IIuoIIh)- (3.17) 


sup 

n>0 




dt 


L°°(0,T;V;) 


This means that the sequence is bounded in L^(0,T; V) n 1T^’^(0, T; V') which is com¬ 
pactly embedded in L^(0,r;H) (Theorem 2.1, page 271 from [23]). Hence, there exists a 
subsequence w^/(w) that converges strongly in L^(0,T;H) to some u‘^{uj) which is also a weak 
limit in L^(0, T; V) and a weak* limit in L°°(0, T; H) and G L^(0, T; V) n L°°{0,T; H). 
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We also have from (3.14) that 


zUu,t) = -l 
1 /■* 

will converge to 2^(0;, t) = - / 

^ Jo 


£ JO 

in C([0,T];H). 


We now pass to the limit when n' —>■ 00 in the system (3.14) pointwise in w G 14. We 

/ 

integrate the first equation over [0,rl, use the convergences of the sequences uL and , so 

ot 

we get that for every (j) G V„ 


lim 



dui, 


and 


TV 

n'—>00 ./q JD 


f dvF 

{s, x)4>{x)dxds = / / -^{s, x)4>{x)dxds 

Jo Jd 


lim 

n'—>00 



Vu^,{s)V(j)dxds = 


0 Jd 



Vu^{s)'V(j)dxds. 


0 Jd 


Also 


a^{z^, + w^)u^^i(j)dxds — 


0 JD 


a'^{z^ + w^)u^4>dx 


0 JD 


a‘^{zl^i + w^){u'^, — u^)4>dxds 


0 JD 


+ 


{ct‘^{z^i + w^) — a^{z^ + w^)) u’^J'dxds 


C 


10 JD 
,£^2 


< 


< 


(u^/ — u^fdxds + C 


/o JD 
, e ||2 


\z^, — z^\ \u^\\(t)\dxds < 


0 JD 


yJ\\i,2(^Q^rp.-^-^ + C j ||z^/ ^;‘'||h||'w''||h||'^||l°o(d)3- 


(3.18) 


so we obtain that 


lim / / {z^, + w^)u^/4>dxds = / / a^{z^ + w^)u^(j)dxds. 

n '^00 Jq Jd Jo Jd 

We use the convergence for and obtain in the limit: 

n -^{t)(l>dxds + f f Vu^{t)V(l>dxds + f f {z^ + w^)u^(pdxds = f f f{t)(l)dxds, 

) Jo Jd Jo Jd Jo Jd 

dz^ 

u^(0) 

I ^^0) 


= u\ 


0> 


= 0 , 


(3.19) 

for every (j) G V„, so by density it is true for any 1 ^ G V. Now, let := z^ + , then we 

deduce that (u'^,u^) is a solution for our initial system in the sense given by (3.2) and (3.3). 
The solution (u^,v^) is measurable as the limit of the Galerkin approximation (n^,u^) which 
is measurable by construction. Furthermore, given the uniform estimates for Uq it is easy to 
obtain from (3.15)-(3.17) the estimates (3.4)"(3.6). 

Now, we prove the uniqueness. Let us assume the we have two solutions {uf, uf} and {^2,^1} 
for the system. Then, 

[ — u\{t))(j)dx + f f {’Vu 2 {s) — Vul{s))'V4>dxds= f f {a^{vf)ul — a^{v 2 )ul)(l>dxds, 

Jd Jo Jd Jo Jd 
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and 

vm - vl{t) = - [\uUs) - uUs))e-^^-^'^/^ds. 
£ Jo 

we take cj) = u^ — and we get: 


Ij 


{u%{t) — u\{t))‘^dx + 



lo Jd Jo Jd 


C 




(Vm|(s) — V uKs))"^ dxds = 

0 Jd 



0 Jd 


|n| — — u\\dxds < 


C 


C / Iki “ ^lllL2(D)3||ni||i4(£))3||tt| - tti||x,4(^)3 < 

Jo 

\ 1/2 / /.t ^ ^ 1/2 


I I|ll2 11 iIIl 2(£))3 ||'ltl||^4(£))3dsj 


'0 


ll^l 11 iIIl4(£))3 ] 


< 


rt 


C / ||n" -""2 

Jo 


^iIIl2(D)3|I^11iIIl2(D)3x 31^S + - / ||Vn| — Vttf ||^2(£))3 x3 

•' 0 


where we used Holder’s inequality and the imbedding of H^{Df into L^{D)'^. Now, 

1 

||u|(t) - u^(t)||L2(^)3 < - ||u|(s) - uf(s)||i2(^)3e"(‘"*)/^ < T sup ||u|(s) - u5(s)||i2(^)3, 

£ Jo se[o,t] 


so we obtain: 


sup ||ll|(t) -<(/)! I < 

s&[0,t] ' 


CT / sup ||ii|(r)-u5(r)||i2p)3||Vuf(s)||^2(^)3x3ds. 

Jo rg[ 0 ,s] 


We use Gronwall’s lemma for the function sup^gjg ^] ||u 2 (^) “ ^i(^)lli (£))3 obtain that: 


CT ||Vuf (s)|||2(£,)3x3ds 

sup ||u|(t) - < ||n|(0) - Ui(0)||2(^pe Jo 

se[o,t] 

which gives the uniqueness and this completes the proof. 


□ 


Theorem 3.2. Assume that the initial conditions Uq are uniformly bounded in V. Then the 
solution will satisfy the improved estimates: 


sup l|ll^llL°°m;L2(o,T;H'2(£))3)) < Ct, 
e>0 

sup ||u^||Loo(n;C([o,r];V)) < Ct, 
£>0 


and 


sup 

£>0 


du^ 


dt 


< Ct- 


L°°(n;i2(o^r;H) 


(3.20) 

(3.21) 

(3.22) 


Proof. To show these estimates we go back to the Galerkin approximation used to show the 

_ Qyf' 

existence. In the system (3.14) we take (f = ~^i^) und get 


'D 


dt ^ ’ 


f dip 

dx+ V<(t)V^(t)dx < G 
Jd ot 


dt ^ ^ 


L2{Dp 


(|I/(^)IIl2(D)3 + I|u^(^)IIl2(D)3) ■ 
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We integrate on [0, t] and use the estimates already obtained for to get: 


Out, 


dt 


is) 


ds+||V<(t)||^<||V<(0)||^ + C / 

H Jo 


dt 


H 


and from here 


and 


sup sup / 

£>0 n >0 Jo 




dt 


ds < Ct, 


H 


supsup sup ||VM^(t)||H < C't, 
£>0 n>ote[o,r] 


which will give us by passing to the limit on the subsequence uL (3.22) and 


sup ||'W*^||l°°(D;L°°(0,T;V)) < C't, 
£>0 


We use now the first equation from (1.2) and the regularity theorem for the stationary Stokes 
equation from [23] to obtain (3.20). We get (3.21) by using Lemma 1.2, section 1.4 from 

□ 


4. The fast motion equation 


In this section, we present some facts for the invariant measure associated with (4.1) and 


(4.1) 


introduce an auxiliary eigenvalue problem. We consider the following problem for fixed ^ G 

dv^ =— ^)dt + ^/QdW, 
i;(0) = r]. 

This equation admits a unique mild solution v^{t) G L^(n; C{0,T; L‘^{D)^)) given by: 

y^t) = r]e-^ + ^{1 - e-^) + [ y/QdW. (4.2) 

Jo 

When needed to specify the dependence with respect to the initial condition the solution will 
be denoted by v^’^{t). The following estimate can be derived for u^’^(t). 


Lemma 4.1. 


E||u^’^(t)|||2(£))3 < ||??||i2(£))3e * + ||^||i2(£))3 + TrQ. 
Proof. It’s enough to use the Ito formula for I|u^’’^(^)lli2(£))3- 


(4.3) 

□ 


4.1. The asymptotic behavior of the fast motion equation. Let us define the transition 
semigroup associated to the equation (4.1) 

P/T(r/) = E^{v^’^{t)), (4.4) 

for every T G and every rj G L‘^{D)^. It is easy to verify that P^ is a Feller 

semigroup because P a.s. 


< e - r]2\\L2(Dp. 


(4.5) 


We also denote by the associated invariant measure on L‘^{DY . We recall that it is invariant 
for the semigroup P/ if 

[ P^^{z)d^^^{z) = [ ^{z)d^iH^), 

Jl'2(dP Jl'2{d)2 
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for every G It is obvious that is a stationary gaussian process. The equation 

(4.1) admits a unique ergodic invariant measure that is strongly mixing and gaussian with 
mean ^ and covariance operator Q. All these results can be found in |10] or [6]. 


by 


We denote by the Kolmogorov operator associated to the semigroup P/, which is given 


f 9T 1 

= —{ri)- [-7] + 0 + 


Q 


drp 


for every 'h G C'|(L^(I1)^). 

For A > 0, let us introduce the eigenvalue problem associated to the operator C^\ 

AT(r/) - = F(r/), 


(4.6) 


(4.7) 


for 7] G LP‘{DY and F G C'^(L^(iA)^). The equation (4.7) has a unique strict solution (see 
Section 9 from M)- Moreover the solution is given by (see Theorem 3.8 in [8]): 


^iv,0= I e-^^P,^F{r,)dt. 


5. Passage to the limit 


The main goal of this section is to pass to the limit in the system (1.2) when e 
We introduce the following averaged operators: 

FF ■. L^{Df ^ [ a^(r/)d/r«(7/) 


0 . 


«(0 = / 
Jl 


L2(D)3 


ot{y,z)dy ) d^^{z). 


’Y 


(5.1) 


(5.2) 


We remark that as an operator from L^{D)^ to L?‘{D)^ is Lipschitz and is sepa¬ 

rable, so Pettis Theorem implies that : Lp‘{D)^ —)• L?‘{D)^ is measurable. The boundedness 
of implies the integrability with respect to the probability measure so is well dehned 
(see Chapter 5, Sections 4 and 5 from [28] for details). The same considerations hold also for 

the operators a e L\Df ^ a(y, e L-(Df, so S is also well defraed. 

Our main result is given by the next theorem. 

Theorem 5.1. Assume the sequence Uq is uniformly bounded in V and strongly convergent 
in H to some function uq, and Vq is uniformly bounded in L^(r2; L^(iA)^). Then, there exists 
u G L^(0,T; V) such that converges in probability to u in L^(0, T; V) and u is the solution 
of the following deterministic equation: 


du 

= Art — Vp — a{u)u + / 

in 

D, 

divlt 

= 0 

in 

D, 

u 

= 0 

on 

dD, 

u(0) 

= uo 

in 

D. 


(5.3) 


Let us explain the main ideas involved in the proof of this convergence. The uniform bounds 
for provided by Theorem 3.2 imply that the sequence is tight in L‘^{0,T; Hq{D)^), so there 
exists a limit u in distribution. We apply after that Skorokhod theorem to get another sequence 
dehned on some probability space 11', with same distribution as that converges for a.e. 
w G O' to some u in L^(0, T; 77^(71)^). We show that u is deterministic and get an equation 
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for it by passing to the limit in expected value in the variational formulation. More precisely, 
we prove first that: 

rT r 


lim E 

S —^0 


= 0 , 


/ / {v^{t))u^{t) — a{u{t))u{t)) (j)il){t)dxdt 

/o Jd 

for every (/> G V and any i/: G (^[O, T] with = 0. We rewrite it as: 

fT /■ 

(a^{v^{t))u^(t) — a{u{t))u(t)) (f>xjj{t)dxdt = 5f + 5 | + 5|, 


(5.4) 


10 Jd 


where 


= 


JO 

rT 


5| = 


f (a^(v^(t)) — a^{u^{t)))u‘^{t)4>'il>{t)dxdt, 
JD 

/ {a^{u^{t))u^{t) — {u{t))u{t)) (j)'tp{t)dxdt, 


and 


10 JD 
rT 


= f f {a‘^{u{t))u{t) — a{u{t))u{t)) 4>il}{t)dxdt. 

Jo Jd 


Thi s convergence requires several preliminary steps. The first step is performed in Subsection 
where we get uniform estimates for and ( 77 , g), where T*" : Lr{Dy x 


5.2 
W{Df 


solves 


dr] 




c(e) 4 '*'(r 7,0 = / (a^r]) - a‘^{^)) ^cpdx. 


(5.5) 


ID 


for c(e) > 0 to be chosen later. This is an eigenvalue problem of the type (4.7) so the solution 
is given by 


4/^ 


poo p 

{V,€)= i / {a^{-)-a^{i))i(t>dx{ri) 

Jo IJD 


dt. 


(5.6) 


The second step is done in Subsection 5.3 where we prove the convergence to 0 for Sf by 
applying Ito’s formula to 4'^(u^(t), tt^(t)) and getting an expression for Sf in terms of 41^ and 
its first order Frechet derivatives. The result is contained in Lemma 15.61 In Subsection 15.41 we 


do the third step, the convergence to 0 of S^, which is showed in Lemma 5.7 


The sequence given by Skorokhod theorem converges a.s. to u strongly in L?‘{Q,T] Hq{D)^) 


so 


lim 
£—^0 


J j (u^{t) — u{t)^ cjrij'(t)dxdt — J j [ Vu^ — Vu ) 'V(f)'ilr{t)dxdt 


= 0, a.s. (5.7) 


The equations (5.4) and (5.7) imply that u satisfies almost surely the variational formulation 


associated with (5.3), so u and u are deterministic and as a consequence the convergence of 


the sequence to u will be in probability. Before proceeding with the proof of Theorem 5.2 


let us first study system (5.3) 


5.1. Well-possedness for the averaged equation (5.3) 


Theorem 5.2. Assume f G L^(0, T; L^(Z9)^) and a G C^(M^) is a positive function. Then, 
for any uq G H the system ( |5.3[ ) admits a unique solution u G C'([0,T];H) n L^(0, T; V) with 
dn 

— G L2(0,r;V') in the following sense: 


dt 


/ u{t)(f>dx — / UQcfdx + / / Vu{s)Vfdxds + / / a{u)u4>dxds = / / f{s)(j)dxds, 

Id Jd Jo Jd Jo Jd Jo Jd 

(5.8) 























HOMOGENIZATION OF BRINKMAN FLOWS IN HETEROGENEOUS DYNAMIC MEDIA 


13 


for every t £ [0, T] and every ^ G V. Moreover, if the initial condition uq G V, then u has the 

dn 

improved regularity, u G L‘^{0,T; n C'([0, T]; V) and — G -L^(0, T; H). 


Proof. The proof of existence of solutions is similar to the proof of system (1.2), using a 
Galerkin approximation procedure. The finite dimensional approximation Un will exist as in 
Theorem 13.11 and will solve 

f dun, r r r 


. '-{t)(j)dx+ / V Un{t)'V (pdx + / a{Un)Un4>dx = / f{t)(j)dx, 

Jd Pi Jd Jd Jd 

for every f G (^([O,T], V„), and Un(0) = li-nUo- We take (p = Un{t), and get: 


(5.9) 


'D 


du„ 

1h 


{t)Un{t)dx+ / \\VUn{t)\\'^dx < / f{t)Un{t)dx 


ID 


ID 


d 


gJ|'Un(t)||j^2(£))3 < II/(^)IIl 2(£,)3 + ||^n(t)||/,2(£))3 
Wn{t)\\i^2(^D-j3 < c + f \\Un{s)\\j^2(^£,'j3ds. 


We use Gronwall’s lemma and get: 


and from here we also obtain 


and 


sup ||n„||c([o,r];H < Ct, 
n>0 


sup ||VM„||x,2(o_'r;L2(D)3x3) < Ct, 
n>0 


sup 

n>0 


dUn 


dt 


< Ct- 


L2(0,T;V') 


(5.10) 

(5.11) 

(5.12) 


So there exists a subsequence Un' and a function u G L°°(0, T;H) n L^(0, T; V) such that Un' 
converges weakly star in L°°(0,T;H) and weakly to L^(0,T; V) to u and also 


du 


dUn' 

1h 


converges 


to — weakly in L^(0,r;V'). We apply again now Theorem 2.1, page 271 and Lemma 1.2 
at 

page 260 from |23j to obtain that Un' converges strongly in L^(0,T;H) and in C([0,T];H) to 
u. We then pass to the limit and obtain that u is a weak solution for (5.3). 

Now, to show uniqueness we assume to have two solutions ui and U 2 in C'([0,T];H) n 
L^(0, T; V) and substract the variational formulations. We get: 


Id 


{U 2 {t) - Ul{t))(l)dx + 



0 JD 


(Vm 2 ('S) — V ui{s))V (pdxds = 



0 JD 


{a{ui)ui — a{u2)u2)(pdxds. 


We take p = U 2 — ui and write 

{a{ui)ui — a{u2)u2){u2 — ui) = —a{ui){u2 — ui)^ + M2(a(u2) — a{ui)){u2 — ui) 

< CU2{U2 — Ul)^. 

We get 

||u 2 (t) - Tl(t)||^2(£,)3 + y ||Vu 2 (s) - Vui(s)||^2(£,)3x3 < c / [ U2{s){u2{s)-Ui{s))‘^dxds 

rt 


0 JD 


< C||u2(s)||i2(o)3||u2(s) - Ul(s)||^4(£,)3ds. 
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after using Holder’s inequality. Interpolation inequality and Young’s inequality give 

||U2(S) -Ul(s)||i4(o)3 < ||W2(S) -Ul(s)||J^i^(^p||VU2(s) - Vui(s)||5^'^2ppx3 

< c(e)||M 2 (s) - Wi(s)|li 2 p )3 +e||Vu 2 (s) - Vui(s)||| 2 (^) 3 x 3 , 
so we obtain for a convenient choice of e 


||u 2 (t) -Ui(t)||^2p)3 < c(e) ||u2(s) - ni(s)||^2(£,)3. 

We get uniqueness from here by applying Gronwall’s lemma. 


Let us now assume that the initial condition uq G V. We use the equation (5.9) with 

dvjTi 


ID 


dUn 

dt 


,du 


dUr, 


. dUr, 


(t) dx+ Vun{t)V^—{t)dx + / a{un{t))un{t)^-{t)dx = / f{t)^-{t)dx, 


ID 


dt 


ID 


dt 


ID 


dt 


we integrate it over [0,r], and use Holder’s inequality: 


(5.13) 


dUr. 


dt 


L2([0,r],L2(D)3) 


+ l|Vttn(Y)||^2(£))3x3 — ||Vtt„,(0)||^2(£))3x3 < C 


dUrt 


dt 


L2([0,T],L2(D)3) 


which will imply that 


dun 

dt 


G L^([0,r],H) uniformly bounded. 


This gives from that Au„ G 


L2([0,T];H) which based on the regularity theorem from the stationary Stokes equation 
implies that Un G L^([0, T], n H) and informally bounded. We deduce from here 

that u G L^([0,r];H) n L^([0, T], and based on Lemma 1.2 page 260 from [23] that 

uGC([0,r];V). 


□ 


5.2. Estimates for We recall that according to (|5.6|) 

\j/'^ 


poo p 

{rj,0= {a%-)-^{0)^4>dx 

Jo UD 


{r])dt, 


defined for a fixed (/) G V n 

Lemma 5.3. 

< c||^||l2||(()||loo (1 + ||r/||i2 + 11^11^2). 


Proof. 


Jo 


(q;*^(-) — a^(0) ^4>dx 


ID 


in) 


(5.14) 


(5.15) 


We denote by F'^{r],f^)= / {r])£^(j)dx, which due to the Lipschitz condition for the function a 

Jd 

is a Lipschitz function of rj with the constant c||^|| 2 , 2(£))3 ||(^||j;^oo (£))3 and Lipschitz with respect 
to f with the constant c|| 0 || 2 ,oo(£)) 3 . 

Using that the measure is invariant, we have that 
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P«F^(r/,e) - [ F%w,C)dfj,Hw) = [ \pfF%rj,C) - PiF^{w,d dfi^w) 

dfi^w) 


JL2(D)3 




using (4.5) 


<ce ^U\\l4Hl°° [ \\{r] - w)\\L 2 dn^{w) 

< ce"*lkllL2||</'llL°° ( M\l^ + [ lhllL2d//^(u;) ) . 

V dL^nr J 

(5.16) 


Again using the invariance of and Lemma 4.1 


/ \\w\\i2d^^{w) = / P^lMli^dn^i w) 


L^{D)3 


= / E\\v^’-{t)h2df,Hw) 

< I c(l + 11 ^ 11^2 + e"*||n^( 0 ,ii;)||i 2 )(i/i^(u;) 

<c(i + ||CIIl2 ) 


(5.17) 


□ 


Now, we estimate the Frechet partial derivatives of 'L'" with respect to rj and 

Lemma 5.4. 




dr] 


iv,0 


L‘2{D)3 


< c\\C\\L^{D)3 Ill’llL°°{D)^- 


(5.18) 


Proof. Let us show first that for every t > 0, 




f da^ 


\ dr] 


which means that 


(5.19) 


lim 

Il'“'lli2(£))3^0 


P«F-(-, 0(r/ + w)- PiF%;0iv) - e-*E (n^’"(t))) ^ • </>) wdx 


\M\l'^(D)^ 
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We have: 


+ w)- 0(1/) - ^ e-^E e • <)>) wdx 

0 - e-*E e • </>) wdx 


=E F 


=E 0 - F%v^’^{t),^) - e“‘ y ^ ■ wC ■ cjidx 

j (t)^ — W • e~^w^ ^ ■ (pdx'^ . 


But from (4.2), = v^''^{t) + e so 

Pid^'i-,0iv + w) - PfF%-,^){i]) - e"*E ^ ■ (jidx 

J (^v^’'^{t) + e~^v?j — W 


e w O 4’dx 


Assume that (5.19) does not hold, so there exists a sequence Wn G L?‘{D)^ with the norms 
converging to 0, such that 


lim E 


q,£ (t)€>’?(t) + e~^Wn) — cF (i;^’^(t)) —— (i;^’^(t)) • e~^Wn ) C ' 4>dx 


ID 


\w. 


i|Il2(D)3 


^0. 


We can assume w.r.g. that Wn —?■ 0 pointwise in D. This implies, by the differentiability of 
the function a that 

rlni^ \ 

qS _(_ e~^Wn) — oF {v^’^{t)) - — ■ e~^WnJ C ■ 4> 


\w. 


i|Il2(D)3 


converges to 0 pointwise in D. We get a contradiction after applying the dominated conver¬ 
gence theorem. 


From (5.19) we obtain by differentiating inside the integral that 


drj 




(5.20) 


which implies (5.18[). 


Lemma 5.5. 




di 


iv,0 


LpD) 


3 ^ c(e) 


(l + ll^llL2(i5)3) ||0||x,oo(o)3. 


Proof. Similarly as in the previous lemma we show that 

Ep5F'(.,{)(„) = E ((1 - e-‘) ^ (-41 +a’ 


□ 


(5.21) 


(5.22) 


and 


d 


da^ 


op, F%w,£,)dfj,^{w)= a%w)(j)dn^{w)+ {w) {w). (5.23) 

Jp{DF JP{Df Jp{DF 
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We write: 

+ z) - + z) - P^F%;C + z) + + z) - 

= E e + z)- F^v^’^it),^ + z) + F^v^’^it),^ + z) - F%v^^^{t),o) 

We have: 

P/+^F^(-, e + ^) - p!f%;C + z)- 1^(1 - e-*)E e • </)) Zdx 

=E ^ + z)- F%v^’^{t),C + z)') - ^(1 - e'^E ^ ' </>) ^dx 


=E ( F%v^+^'^{t),^ + z)- F%v^’\t),^ + z)-il- e-*) 


Id 


da'^ 


(^v^’^it)^ ■ z^ ■ (f>dx'^ 


da^ 

But from (4.2), + (1 — e~^)z, so 


=E { I [a^ (— W 




(1 — e *) 2 : K • (pdx 


P«+"F-(-, e + ^) - Pt^F%;^ + z)-il- e-*)E 
=E 




Id ■z^-cpdx'^ 

j (v^’^^it) + (1 - e~^)z'^ - W (v^’'^it)^ - ■ (1 - e~^)z ) ^ • 0dx ) . 

We show that 


limE 

z^O 


da^ 

q,£ + (1 — e~*)z) — (u^’'^(t)) —— {v^’^{t)) ■ (1 — e~'^)zn ) ^ • (j)dx 


'D 


z 


IIl2(_D)3 


= 0 . 


(5.24) 


By contradiction assume that there exists a sequence Zn € L'^{D)^ with the norms converging 
to 0 , such that 


lim E 


Bcy^ 

a'^ + (1 - e~^)zn) - W (u^’^(t)) - ■ (1 - e~^)zn ) ^ • Hx 


/O. 


a 


Jd lkn||L2(£))3 

We can assume, by passing to a subsequence, that Zn converges also pointwise to 0. But by 
the differentiability of the function a, we have that 

+ (1 - e~^)zn) - • (1 - e~^)zn] C ' 4> 

\\^n\\p{D)3' 

converges to 0 pointwise in D. We get a contradiction after applying the dominated conver¬ 
gence theorem. 

We use similar arguments and calculations to prove that 


lim 

2:—>0 


zdx 


P!F%;C + z) - P!F%;C) - 

lkllL2(i5)3 

The limits (5.24) and (5.25) imply (5.22) which implies that 

Pt P^i'^Oiv) < c||(/>||£,oo(£,)3(||^||i2(£))3 + 1). 


= 0 , 


Ap? 


(5.25) 


(5.26) 


L2(P)3 

















18 


H. BESSAIH, Y. EFENDIEV, AND F. MARIS 


□ 


Now we compute — [ a^{^)^(l)dx = 7 ^ / [ {z)dfi^{z)^(j)dx and we show: 

c'? Jd <^4 Jd Jl‘2(d)^ 


d f — 




da‘^ 


'D 


a^{^)^4>dx= I '^^^{z)d^^{z)^(j) + [ a^{z)d^^{z)(j). 

L2(D)3 or] Jl^(D)3 


(5.27) 


The measure ]j,^ having the mean ^ and covariance Q satisfies ]r^{B) = ]i^{B + for every 
borel set B G Lp‘{D)^. So using the change of measure: 

^ [ [ a^{z)dn^{z)^(pdx = ^ , , 

01; Jd Jl^{d)^ 01; J^, Jl2{d)^ 

Now (5.27) follows using the same technique as for (5.24[) and (|5.25[) and we also get: 


a^{z)dfi^{z)^4>dx = — I I a^{z + ^)d]/'{z)^(pdx. (5.28) 




< c\\4>\\l°°{D)3{\\(\\l^{d)3 + !)• 


(5.29) 


We use the estimates (5.26) and (5.29) in the expression for and obtain (5.21). 


5.3. Convergence of Sf. We start by applying Ito’s formula to (v^(t), (t)), for a fixed 

^ E VnL“(T))3. 

(u^, u^) + ^iv%s),u%s)) • {Au%s) - a%v%s))u%s) + fis))ds 

+ r + y^%s))ds ^(y^(s),U^{s)) • ^dW{s) 


52^6 


JO 

1 ^ ^ 

+ — Tr Q ^ ^ (v^(s),u‘^(s)) ds 
Jo ^ 

= + I* ^(u-(s),u^(5)) • iAu%s) - a%v%s)Xis) + fis))ds 

1 /"* 

+ - / £“'WT"(u"(s),u"(s))ds 
£ Jo 

1 7* 


(5.30) 


+ ^ / {v^{s),u%s))-yJQdWis). 

Jo or] 

We use the expression for the Kolmogorov operator to get 

^%v%t),u%t)) = + [ ^^{v%s),u%s)) ■ {Au%s) - a%v%s))u%s) + f{s))ds 

Jo Ol; 

_l_ f ^^[v^[s),u^{s))ds + - f f {a^{v^{s))ds — a^{u^{s))) u^{s)(j)dxds 

^ Jo ^ Jo Jd 

Ve Jo or] 

f f {a^{v^{s)) — a‘^{u^{s)))u^{s)(l)dxds = £'i!^{v^{t),u^{t)) — £^^{vq,Uq) 

Jo Jd 

- ^ f • [Au^{s) - a^{v^{sW{s) + f{s))ds (5.31) 

Jo Oi; 

-c{£)[ ^^v%s),u^{s))ds - ^/£ [ ^^{v^{s),u^{s)) ■ ^/QdW{s). 

Jo Jo Or] 
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Lemma 5.6. Let us define by Sf the integral 



{a^{v^{t)) — a^(u^(t))) {t)^^'{t)dxdt. 


0 JD 


Then for a particular choice of the sequence c(e) we have the following convergence: 

lim E IS'f I =0. 

£—>0 


(5.32) 


Proof. From the equation (5.31) we get: 

rt 


{a'^{v^{s)) — a^{u^{s))) u^{s)(pdxds 

/o JD 


+c(e) / \^^{v^{s),u‘^{s))\ds 

Jo 




dr] 


(u''(s),u^(s)) • y/QdW{s) 


(5.33) 


+£ 


10 




{v^(s), (s)) ■ (Au^(s) — a^(u^(s))M^(s) + f{s)) 


ds. 


We use for the function the estimates (5.14), (5.18) and (5.21) and get: 



(a^(v^(s)) — a^(u^(s))) u^(s)fidxds 


0 JD 


< 


ce||u'^(t)||i2p)3||(/)||ioo(£))3 (l + \\v^{t)\\L‘2(Dfi) + ce||uollL2(D)3 II</'IIl°°(d) 3 (l + II^oIIl2(d)3) ds 


lo dv 

rt 


{v^{s),u%s)) ■ y^dW{s) 


+c(e) f ce||M^(s)||i2(£))3 ||(/)||x,oop)3 (l + ||u^(s)||x,2(£))3) ds 
Jo 

f 

Jo 


(u^(s),u^(s)) 


•||(Au^(s)-W(u^(s))u^(s) + /(s))|L2(^)3 


-\-S ^ ,1 ' 

.L2(D)3 

<C£\\4>\\l^]DP sup ||w^(i)||L2(£))3(l + ||l^^(i)||L2(£))3 + 11 Uq 11^2(15)3) 
te[o,r] 

l'^^{v%s),u%s))-^dWis) 

+cec(e)||(/>||£,oop)3 sup ||ii^(t)||L2(^)3 / (l + ||u^(s)||£,2p)3) ds 
te[o,T] Jo 


-\-s 


[q c{^ WJ’Wl^IDP ■ (ll(^^(s)||_H-2(£))3 + II/(s))IIl2P)3^ 


(5.34) 


But, Ito’s isometry gives us that 

ft 


E 


to 


dg 


(u^(s),u^(s)) • y/QdW{s) 


< c 


< 




dg 


{v%s),u%s)) 


\ 1/2 


L2(D)3 


J 


^^=II</’IIl“(d) 3 (l + I|ii*^('S)||i, 2(£))3) ds 


1/2 
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So, given the estimates for u^, we obtain that 

E|5f| <CTeE(l+ sup \\v^{t)\\L^^D)^) + + Ct^. 

te[o,T] V 

Now using the estimate for and choosing c(e) = y/e we get (5.35|). 


□ 


5.4. Convergence of S^. 


Lemma 5.7. For fixed u G L°°(0; (^([O, T]; L^(Z))^)) and fi G let us define by the 

integral / / {a^{u{s))u{t) — a{u{t))u{t)) (l>'fi{t)dxdt. Then: 

Jo Jd 


lim E IS*!! =0. 
£—>■0 


(5.35) 


Proof. For any t G [0, T] consider the sequence of functions Ff : Lp‘{D)^ —^ 

Ft{^){x) = {^,z{x)^ - a {y,z{x))^ u{t,x). 

We show now that for any z G LP‘{DY, for every t G [0,r] and a.e. u £ 0., Ff{z) converges 
in L^{D)^ to 0. We we fix u and t and let Zn and Wn two sequences of continuous functions 
converging in L^{D)^ to z and u{t). We use Lemma 1.3 from [T] and obtain that the sequence 

Ff{x) = ^ —, Zn{x'^ — J a{y, Zn{x))^ Wn{x) converges when e — )> 0 to 0 in L‘^{D)^. 

But 

\Ff(x) - Ff{z){x)\ < c\Wnix) - u{t,x)\ + c\Znix) - z{x)\, 
based on the Lipschitz condition and boundedness for a. We deduce that Ff{z) converges in 
LfilyDY to 0. The sequence being also uniformly bounded by ||h||j;^oo(Q.(^([o,r];L2(D)3))) Vitali’s 
convergence theorem implies that the sequence of the integrals with respect to the probability 
measure on also converge to 0 in L\Df-. 

lim/ Ff{z)dfPJ^dz = 0in 

which can be rewritten as 

lim a^{u(t))u(t) — a(u(t))u{t) = 0 in LfiiD)^. 

£->0 

This implies that E a.s. and for every t G [0, T] 

lim / {a‘^(u{t))u(t) — a(u(t))u{t)) fiif'iPdx = 0, 

with the sequence being also uniformly bounded by c||1Z||x,oo(q ||<('||lo°(£))3 ||V'1Il°°[o,t]- 

We apply the bounded convergence theorem and integrate over x [0, T] to get the result. □ 


5.5. Proof of Theorem 15.IL 


Proof. The uniform bounds (3.20), (3.21) and ( 3.22[ ) hold for u^. Then, for a.e. a; G 12 the 
sequence is included in a compact set /C of L^(0,T; V) and the sequence is tight. Then, there 
exists a subsequence and a random element u G L^(0,T;V) such that converges in 
distribution to u. Skorokhod theorem gives us the existence of another subsequence and a 


sequence u^" with the same distribution defined on another probability space 12 that converges 
pointwise to some TL, a random element of L^(0, T; V) with the same distribution as u. It follows 
from here that u belongs to /C a.s. so n G L°°(12, L^(0, T; V)) and u G L°°(12, L^(0, T; V)). We 
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use the variational formulation ( |3.2[ ) for u^'' with a test function i;/) G V n multiply it 

with Tp', where G with V’(T) = 0 to get: 


f f {t)dxdt — [ f 

/o Jd Jo Jd 

rT 


10 


Uq 4>'ip {t)dxdt — / / {t)VJ)'p{t)dxdt 

^ (5.36) 

J {v^ {t)(j)'tp{t)dxdt = — f [ f{t)4>'ip{t)dxdt. 


10 Jd 


We now show that 


lim E 

£"-S>0 


/7 

/o Jd 


W {v^ it))u^ (t) — a{u{t))u{t)^ (j)'tp(t)dxdt 


= 0 . 


(5.37) 


We write: 


W'V'V"" - Hu)u = a^"{v^"X" - a""(n"'>"" + a""(n"")u"" - W"(u)n"" 
+ a^" iu)u^ — a.^'' (u)u + W” (u)u — a{u)u, 

r [ a^"(v^"it)X''it)(j)i;{t)dxdt = sf + sf + sf, 

Jo Jd 




SO 


where 


SI = / {a‘^{v'^{t)) — a^{u^{t)))u^{t)(f)'ip{t)dxdt, 

Jo Jd 

f f {a^{u^{t))u^{t) — a^{u{t))u{t)) (j)'ip{t)dxdt^ 

0 Jd 


Sl = 


and 


5f = 


(5.38) 

(5.39) 

(5.40) 


/ / {a'^{u{t))u{t) — a{u{t))u{t)) (j)'4i{t)dxdt. 

lo Jd 

Given the same distribution for the sequences and u^, for u and u and the uniform bounds 
in L^(0, T; V) we have: 

= E / f (a^"{u^''{t))u^"{t) — a^"{u{t))u{t'^ (l)'p{t)dxdt 


E 




< E 


+ E 


Jo Jd 

/7 

JO Jd 

fl 

lo Jd 


(u^"{t))u‘^"(t) — a^"(u'^”{t))u{t )) (f)'ip{t)dxdt 


a^"{u^"{t))u{t) — a'^”{u{t))u{t )) 4>ip{t)dxdt 


The function a*" is Lipschitz, with the same constant as a and also bounded so we get: 


E 




< c||(/)||ioo(£))3||V'||ioo[0,T]IE 


f" — 

— u 


li(0T;V) 


The pointwise convergence of u^" to u in L^(0,T; V) and the uniform bounds imply that 


lim E 
£"^•0 


S% 


= 0 . 


This, together with the limits given by the Lemmas 5.6 and 5.7 give (5.37). (5.37) and (5.36) 
imply that 

E y J (t) — Uq pip'{t) — V{t)'Vpip{t) + f{t)pip{t) — a{u{t))u{t)pip{t)^dxdt 

Now, using the fact that and are equally distributed implies that 


0. 
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E 


ID 


u^"{t) — Uq (pip'it) — Vu^''{t)V<j)'ip{t) + f{t)(p'tp{t) — a{u{t))u 


dxdt 


The sequence above converges in L^(0) to 0 but also pointwise in Q to 


Id 


u{t)(f)'ip'{t) — uo4>ijj'{t) — Vu{t)V4>'ijj{t) + f(t)4>'tp(t) — a{u{t))u{t)(j)il>{t)J dxdt, 


which means that u i s po intwise the weak solution of the deterministic equation (5.3) which, 
according to Theorem . 


5.2 


0 . 


has a unique solution, so u and u are deterministic. Then, the whole 
sequence converges to u in distribution, and since u is deterministic then the convergence 
is also in probability see m Theorem 18.3. □ 
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